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Abstract. We review the main stellar features observed in starburst spectra from the 
UV to the near-IR and their use as fundamental tools to determine the properties of 
stellar populations from integrated spectra. The origin and dependence of the features 
on stellar properties are discussed, and we summarise existing modeling techniques used 
for quantitative analysis. Recent results from studies based on UV, optical and near- 
IR observations of starbursts and active galaxies are summarised. Finally, we briefly 
discuss combined starburst + photoionisation models including also observations from 
nebular emission lines. The present review is complementary to the recent summary by 
Schaerer (2000) discussing more extensively nebular analysis of starbursts and related 
objects. 



1 Introduction 

The analysis of distinct spectral features in integrated spectra is at the base 
of numerous investigations on the stellar content of distant galaxies. Indeed, in 
addition to the overall continuum spectral shape, stellar absorption (and rarely 
also emission) lines carry crucial information on the presence of stars of various 
spectral types and luminosity class, and allow thus in principle to "decompose" 
the integrated galaxy spectrum in its stellar constituents, and to determine their 
fundamental properties such as ages, IMF, the star formation history etc. 

In objects such as active galaxies, where non stellar emission processes are 
thought to contribute to the emitted light, the study of possible stellar features 
allows one to constrain the relative stellar contribution, and thus to determine 
the efficiency of various emission processes (e.g. stellar versus non-stellar). 

These basic properties illustrate the interest of spectroscopic studies of stellar 
features in starbursts and other galaxies. 

The aim of the present review is to discuss the main stellar features observed 
in starburst spectra over the entire spectral range where such features are de- 
tectable, i.e. from the UV over the optical to the near-IR. At longer wavelength 
dust emission dominates and the stellar continuum and associated lines is not 
detectable anymore. 

An "inventory" of the strongest stellar features is given for each spectral 
domain and recent results in the respective fields are summarised. In the last 
Section, we also briefly discuss the use of combined stellar and nebular emission 
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line (hereafter EL) analysis for the study of stellar populations in the otical and 
IR. 

The current review is complementary to a recent review Schaerer (2000) 
discussing new developments in multi-wavelength modeling tools, the current 
status of ionising fluxes from massive stars, and their importance for EL analysis 
of starbursts and related ojects. 

2 UV features 

The UV spectral range (~ 1000 - 3000 A) of starbursts is rich in stellar lines 
originating in early type stars (mostly OB, also Wolf-Rayet), it contains few or 
weak nebular lines, and rather numerous interstellar (IS) absorption lines. 

A rough inventory of the strongest stellar lines observed mostly in the ~ 1200 
- 1800 A range follows (cf. the detailed work of de Mello et al. 2000, hereafter 
DLH00): 

• Well known stellar wind lines (P-Cygni or EL) from O and Wolf-Rayet (WR) 
stars) are Si iv 1400, C iv 1550, N v 1240, He n 1640, and N IV1720. The 
following synthesis models include at least partly these lines and discuss 
their behaviour: Sekiguchi & Anderson (1987), Mas-Hesse & Kunth (1991), 
Fanelli et al. (e.g. 1987, 1992), Leitherer, Robert et al. (1993-2000). 

• Other wind lines blueward of Lya include O vi+Ly/?+C II 1010-1060, dis- 
cussed and modeled by Gonzalez Delgado et al. (1999), and potentially lines 
of C, N, P, S, and Ar in the range recently observed with FUSE (see Taresch 
et al. 1997, Fullerton et al. 2000). 

• The strongest photospheric lines from OB stars are Si II 1265, 1485, Si iii1295- 
1300, 1417, C ii 1334, 1335, C m 1247, 1427, S v 1501 (sec de Mello et al. 
2000). 

• Other stellar features include "depressions" due to numerous Fe lines (at ~ 
1400, 1600, 1940), and also Fe n 2570-2615 and Mg II 2780-2825 features at 
longer wavelengths (e.g. Robert et al. 1999, Storchi-Bergmann et al. 1995) 

It is important to note that many of these lines can also be formed in the inter- 
stellar medium (cf. Heckman & Leitherer 1997, Sahu 1998). A careful separation 
of the stellar and interstellar component is necessary in many cases (see de Mello 
et al. 2000). The wavelength range between ~ 1800 and 3000 A remains still little 
explored. A similarly detailed understanding of this spectral range is of impor- 
tance for studies of galaxies in the z ~ 1-2 redshift range (see e.g. Campusano 
et al. 2000). 

Given the dependence of the various lines on stellar luminosity, age, and also 
metallicity (cf. below), e.g. well studied for the stellar wind lines (e.g. Walborn 
et al. 1985, Leitherer & Lamers 1991), the features can be used to constrain 
the parameters of the integrated population, such as age, SF history, and IMF, 
by means of evolutionary synthesis techniques (e.g. Mas-Hesse & Kunth 1991, 
Leitherer et al. 1995). The most up-to-date model suited to such analysis is 
Starburst99 (Leitherer et al. 1999, de Mello et al. 2000). 
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These techniques have been extensively applied to the interpretation of UV 
spectra of nearby starbursts (mostly HST spectra), especially by Leithcrcr, Heck- 
man, Gonzalez Delgado and collaborators (some references given below) and by 
Mas-Hesse & Kunth (1991, 1999). Summarised in one sentence (...) the main 
result of these studies is that all the objects contain young bursts ( < 10-20 
Myr) characterised by instantaneous burst or continuous star formation, the 
distinction being often difficult to draw, which are populated with a rather nor- 
mal Salpeter-like IMF with stars up to M up ~ 60 - 100 M Q . In a recent study 
Tremonti et al. (2000, cf. these proceedings) examine the stellar populations in 
the field of NGC 5253 and find a possible indication for a steeper IMF, although 
other explanations (e.g. age effects) are possible. 

The similarity of the spectra of many high redshift galaxies (e.g. Lyman 
break galaxies) with the local starbursts is now well recognised and offers many 
exciting possibilities. For example, from the beautiful spectrum of the lensed z 
~ 2.7 Galaxy 1512-cb58 of Pcttini ct al. (2000) these authors and dc Mello et al. 
(2000) derive a constant star formation, an IMF slope between Salpeter (2.35) 
and ~ 2.8, and find indications for a subsolar metallicity, in agreement with EL 
measurements from Tepliz et al. (2000). 

Obviously it is of great interest to derive/estimate the metallicity (Z) from 
stellar UV lines. Since, for example, the strength of stellar wind lines depend 
on Z this is in principle possible. This can e.g. be done using the correlation of 
the equivalent width of Si iv with metallicity found by Heckman et al. (1998). 
However, a priori, such a correlation should only be valid in a statistical sense, 
since the line strength also depends on age (cf. Lcithcrer et al. 1995). This 
difficulty should be less if the full wind line profile can be analysed. The inclusion 
of spectral libraries of metal-poor stars in evolutionary synthesis models has 
just been completed (e.g. Leitherer et al. 2000, see also Heap 2000). Alternative 
possibilities to derive the chemical composition include the use of IS absorption 
lines (cf. Pettini et al. 2000), or the use of weak stellar features such as various 
Fe blends known to vary with Z (cf. Haser et al. 1998). 

3 Stellar features in the optical 

The inventory of the strongest stellar lines in the optical is as follows: 

• Broad emission lines from Wolf-Rayet stars of various subtypes (WN, WC) 
are detected in some young starbursts: He II 4686 bump, C iv 5808, C in 
5696, possibly also N in 4512, Si in 4565 (see e.g. Schaerer et al. 1999b, 
Guseva et al. 2000). Synthesis models treating these lines include Cervino, 
& Mas-Hesse (1994) and Schaerer & Vacca (1998). 

• H and He absorption lines from OBA stars have been discussed and modeled 
by Diaz (1988), Olofsson (1995), and Gonzalez Delgado et al. (1999). 

• The Ca II triplet - 8498, 8542, 8662 has often been studied (e.g. Terlevich 
et al. 1990ab, Mayya 1997, Garcia- Vargas et al. 1998). Its origin in both late 
type giants and supergiants complicates a priori the analysis. 
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• Other metallic features and molecular bands originating in stars with F types 
and later are Ca n H+K 39XX, CH G band 4284-4318, Mg 1+MgH 5156- 
5196, Na I 5880-5914, various TiO bands > 6200. These are found e.g. in the 
template spectra of Bica & Alloin (1986 and subsequent papers) for clusters 
and in the starburst spectra of Storchi-Bergmann et al. (1995). 

A complete overview of all starburst studies exploiting these features is not 
possible here. I shall instead briefly summarise recent results on starburst (and 
possibly also AGN) studies using Wolf-Rayet features and metallic lines. The 
contribution of Gonzalez Delgado (these proceedings) illustrates the use of H 
and He absorption lines. 

3.1 WR features as a probe of the most massive stars 

Since WR are the descendents of the most massive stars, detections of their 
features provide the best indication of the presence of massive stars (Minitial > 
25-60 M ) and allow to constrain the upper end of the IMF. A catalogue of all 
known galaxies with WR. detections has been compiled by Schaerer et a l. (1999b; 
on the Web at jittp : //webast . ast . obs-mip . f r/people/ schaerer/ ). 



Starburst— AGN connection The detection of stellar features including the 
so-called WR-bump, UV lines and the Ca II triplet in the Seyfert2 galaxy Mrk 
477 by Heckman et al. (1997) has considerably re-vived this subject, before more 
focused on detections of the Ca II triplet (e.g. Terlevich et al. 1990ab). Since then 
other possible WR detections indicating important massive star populations have 
been made in Mrk 1210 (Sey2, Storchi-Bergmann et al. 1998), Mrk 463E and 
Mrk 1 (Sey2, Gonzalez Delgado et al. 2000), TF 1736+1122 (Sey 2, Tran et 
al. 1999), and in three PG QSO (Lipari et al. 2000). Given the strong nebular 
contamination due to He n A4686, other unambiguous massive star features are 
required to fully clearly establish the presence of massive stars in these objects. 
Such attempts, aiming to detect the WR lines of C iv and/or C in which are not 
affected by nebular contamination, have been undertaken by Kunth & Contini 
(1999) with 2-D spectroscopy in several Sey2. 

From the analysis of WR and H+He absorption lines in their sample of 20 
Sey2 galaxies, Gonzalez Delgado et al. (2000; cf. these proceedings) find that the 
blue and near-UV light of half of their objects is dominated by young and/or 
intermediate age stars. A similar result was found by Storchi-Bergmann et al. 
(1998) on a smaller sample. 

"Normal" starbursts (so-called WR galaxies) Studies on WR galaxies 
(mostly BCD, Irr, spirals) are summarised in the reviews of Schaerer (1999ab). 
Including the detections of spectral signatures from both WN and WC stars in 
a fair number of objects covering a large metallicity range, the following overall 
conclusions emerge from the studies of Schaerer et al. (1999a) and Guseva et 
al. (2000). Except possibly at the lowest metallicities a good agreement is found 
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between the observations and the evolutionary synthesis models of Schaerer & 
Vacca (1998). From this comparison one finds clear indications for short bursts 
(At < 2-4 Myr) in objects with subsolar metallicity, an IMF compatible with 
Salpeter, and a large upper mass cut-off of the IMF, in agreement with several 
earlier studies. In addition, the observed WC/WN star ratios provide new con- 
straints for mass loss and mixing scenarios in stellar evolution models (Schaerer 
et al. 1999a). 

We have recently undertaken a first study of metal-rich starbursts (metallic- 
ities up to ~ 2-3 times solar) with the aim of constraining the upper end of the 
IMF in such environments. From the strengths of the observed WR features we 
derive a conservative lower limit of M up > 30-40 M (Schaerer et al. 2000). New 
observations are being obtained to improve the accuracy of this result. Direct 
studies of the stellar content are of prime importance, also to verify the reliability 
of indirect studies based on nebular line analysis (cf. below). *** re f s *** 

3.2 Population synthesis studies using metallic lines 

Early starbursts studies using metallic lines have mostly concentrated on the 
Calcium triplet (see references above). More complete analysis of starbursts and 
AGN spectra using numerous optical stellar features have recently been pre- 
sented by two groups. Both approaches are based on "classical" population syn- 
thesis using either observed stellar templates (Serote Roos et al. 1998, Boisson 
et al. 2000) or cluster templates (Raimann et al. 2000ab). 

The former authors analyse 5000 - 8000 A spectra of 12 starburst, Seyfert 
and LINERs, and use a synthesis technique yielding a mathematically unique 
solution to determine the relative contributions of different stellar populations. 
Regarding e.g. the importance of super metal-rich stars their results differ from 
other work (e.g. Cid Fernandes et al. 1998, Gonzalez Delgado et al. 2000). 

Raimann et al. have analysed average spectra of H n galaxies, starbursts 
and Sey2 taken from the Tcrlcvich et al. (1991) catalogue. Their main result 
regarding H n galaxies is the finding of significant old (< 500 Myr) underlying 
populations, which also modifies emission line diagnostics based on equivalent 
widths. Their conclusions are supported by the study of selected BCD's by Mas- 
Hesse & Kunth (1999) and by comparisons of predicted and observed EL trends 
in several large samples of H n galaxies (Stasihska et al. 2000). 

4 Stellar features in the near IR 

This spectral range is in most cases dominated by features from late type stars 
(G, K, M) corresponding to red supergiants (RSG), asymptotic giant branch 
stars (AGB), or red giants. Typical ages for the appearance of these stars are > 
10 7 , 10 8 , and 10 9 yr respectively. The strongest features in the K and K band 
are 
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• molecular features, such as CO (6,3) 1.6, CO (2,0) 2.29 and many OH lines 
(all wavelengths given in (im here). 

The following papers (incomplete selection) describe their dependence on stellar 
type and luminosity class and/or provide spectral libraries: Kleinmann & Hall 
(1986), Lancon & Rocca-Volmerange (1992), Origlia et al. (1993), Ali et al. 
(1995), Forster-Schreiber (1998). Recently, some authors (e.g Gilbert & Graham, 
these proceedings) have begun to use also theoretical spectra of cool stars for 
population studies. 

Rather than providing a detailed review of the many studies undertaken in 
this area, I will briefly recall some difficulties affecting near-IR studies of stellar 
populations. While obviously the traditional method of population synthesis, 
decomposing the integrated spectrum in various template constituents, can be 
equally applied to any wavelength range, the near-IR properties predicted by 
evolutionary synthesis models are unfavourably affected by uncertainties in post 
main-sequence stellar evolution and the modeling of cool stars. 

Indeed, it is found that all current non-rotating stellar evolution models pre- 
dict an incorrect variation of the relative red/blue supergiant lifetimes with 
metallicity Z (Langer & Maeder 1995) []. In addition the predicted T c g of RSG 
may be too high (Origlia et al. 1999). This implies, e.g. that the CO features pre- 
dicted by evolutionary synthesis models based on these tracks are weaker than 
observed at subsolar Z (Origlia et al. 1999). Also, the predicted strong metallic- 
ity dependence of colors like V-K (see Cervino & Mas-Hesse 1994) is therefore 
incorrect. Improvements are expected from new stellar models allowing for more 
realistic mixing scenarios including rotation (cf. Maeder & Meynet 2000) and 
better understanding of mass loss in these phases. 

Uncertain mass loss scenarii render the prediction of AGB stars, whose in- 
fluence at long wavelength is non negligible for populations with ages > 1 Gyr 
(e.g. Bruzual & Chariot 1993), rather difficult, as best illustrated by Girardi 
& Bcrtelli (1998). For regions with small masses stochastical effects (Lancon & 
Mouhcinc 2000) and fluctuations of the IMF from the finite number of stars 
(Cervino et al. 2000) also lead to an expected dispersion. While these effects are 
obviously of general nature, they are of particular importance for predictions in- 
volving the dominant contribution from stars with very short lived phases (e.g. 
AGB, WR). 

5 Combined stellar and nebular analysis of starbursts 

In various situations (e.g. starbursts with strong optical EL; IR observations - 
A > 4-10 fim where the continuum emission is dominated by dust and stellar 
features are thus completely absent) analysis of nebular emission lines are of 
interest to constrain the stellar population. 

1 Depending on the adopted set, the predicted RSG/BSG may well be correct for 
a certain metallicity (e.g. Geneva models ok for solar metallicity). The predicted 
metallicity variation does, however, not follow the observed trends. 
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However, given the very nature of nebular physics, the EL are not only sensi- 
tive to the ionising spectrum carrying information of the stellar populations, but 
depend also strongly on the nebular geometry and chemical composition. The 
dependence on these additional parameters (essentially the so-called "ionisation 
parameter" U and composition) render EL studies of stellar populations more 
complex and require thus the use of sufficient observational constraints. 

Recent developments and the current state-of-the-art of stellar ionising fluxes 
forming the input to photoionisation models have been reviewed by Schaerer 
(2000) and shall not be repeated here. In the following we briefly summarise the 
main recent studies undertaken in the optical and IR (cf. Schaerer 2000 for more 
details.) 

5.1 Optical studies 

Recent tailored starburst and photoionisation models are presented in the studies 
of Gratia- Vargas et al. (1997: NGC 7714), Luridiana et al. (1999: NGC 2363, 
2000: NGC 5461), Stasihska & Schaerer (1999: I Zw 18), and Gonalcz Dclgado & 
Perez (2000: NGC 604). Although somewhat different in each study, the general 
approach is summarised in Gratia- Vargas et al. (1997). 

Overall one finds that both the stellar and nebular lines give consistent results 
regarding the main properties of the stellar population, such as age, IMF etc. 
At a more detailed level, however, several of these studies encounter significant 
difficulties (e.g. the temperature sensitive ratio [O in] A4363/5007 is underpre- 
dicted), which indicates that some physical are missing in the photoionisation 
models (Stasihska & Schaerer 1999, Luridiana et al. 1999, also Stasihska et al. 
2000). In short, although most observables can be reproduced by the combined 
starburst and photoionisation models — and the tool can thus be used to de- 
rive SB properties from the EL — one has to conclude that for accurate studies 
relying on nebular lines from H n regions (and presumably also more complex ob- 
jects) some additional physical process(es) (possibly shocks, conductive heating 
at X-ray interfaces etc.) must be taken into account (cf. Schaerer 2000). 

5.2 Starburst + photoionisation models in the IR domain 

Analysis of IR observations (mostly from SWS and LWS on ISO) of starbursts 
based on combined SB + photoionisation models are just beginning to appear in 
the literature. In this context it is useful to keep some intrinsic difficulties in mind 
are. Given the nature of objects and the large apertures involved, the integrated 
spectrum generally includes a large variety of regions. This fact, together with the 
complex geometries involved, render a priori the construction of photoionisation 
models difficult. 

Simple models were constructed for case studies of Arp 299 and M82 by 
Satyapal et al. (1998) and Colbert et al. (1999) to interpret their LWS (40-200 
ixm) spectra. Colbert et al. (1999) find that the observed EL spectrum of M82 
is compatible with an instantaneous burst at ages ~ 3-5 Myr, a Salpeter IMF, 
and a high upper mass cut-off. Surprisingly, inspection of models with similar 
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ingredients (cf. Stasiriska & Leitherer 1996), show that the shorter wavelength 
data (see Genzel et al. 1998) is clearly incompatible with the Colbert ct al. 
model predicting too hard a spectrum. In view of the few line ratios originating 
from the H II gas and the large number of free parameters the photoionisation 
model is underconstraincd. A larger wavelength coverage or other constraints 
are required. 

Forster-Schreiber (1998) has described the geometry of clusters and gas 
clouds in M82 by a single "effective" ionisation parameter. This value has been 
adopted as typical for a sample of 27 starbursts in the SB + photoionisation 
models of Thornley et al. (2000). Instead of modeling a simple stellar popula- 
tion their models are based on an ensemble of H II regions following an observed 
luminosity function, which overall leads to a reduction, albeit small, of the hard- 
ness of the ionising spectrum. From the ISO/SWS [Ne m]/[Ne n] line ratios they 
conclude that the observations are compatible with a high upper mass cut-off 
(M up ~ 50-100 M Q ). To reproduce the relatively low average [Ne m]/[Ne n] 
ratio, short timescales of SF are required. More detailed studies including addi- 
tional observational constraints would be very useful to confirm this result. 

A different approach has been taken by Schaerer & Stasiriska (1999), who 
modeled two well studied objects (NGC 5253, II Zw 40) with a fairly well known 
massive star population and existing UV-optical-IR observations. While their 
model successfully reproduces the stellar features and the observed ionisation 
structure of H, He, and O (as revealed from the optical and IR lines), the 
predicted IR fine structure line ratios of [Ne ni]/[Ne il], [Ar in]/[Ar n], and 
[S iv]/[S in] show too high an excitation. The origin of this discrepancy (atomic 
data? separate emission components ? other?) is still unknown. In any case this 
attempt to describe two relatively "simple" objects illustrates the current lim- 
itations and shows that further progress is needed for a proper understanding 
and use of the IR fine structure lines as reliable diagnostics. Improvement is 
expected from multi-wavelength analysis of simpler objects (e.g. Galactic and 
LMC H II regions, PN) and other ongoing work. Such studies should be crucial 
to reliably extend the diagnostic tools to the IR to fully exploit the enormous 
observational capabilities provided by recent and upcoming facilities in probing 
the properties of massive star formation from the local Universe to high reshift. 

Acknowledgements I thank the organisers for this very interesting and stim- 
ulating workshop and for financial support. Part of this work is also supported 
by the INTAS grant 97-0033. 
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